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In this work, the relation between muon production in extensive air showers and features of 
hadronic multiparticle production at low energies is studied. Using CORSIKA, we determine typical 
energies and phase space regions of secondary particles which are important for muon production in 
extensive air showers and confront the results with existing fixed target measurements. Furthermore 
possibilities to measure relevant quantities of hadron production in existing and planned accelerator 
experiments are discussed. 
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1 Introduction 

The energy spectrum and composition of the cosmic rays with energies above 10 15 eV 
are typically derived from measurements of the number of electrons and muons produced 
in extensive air showers (EAS) at ground. However, the results of such a shower analysis 
are strongly dependent on the hadronic interaction models used for simulating reference 
showers (J], Therefore it is important to study in detail the role of hadronic interactions and 
in particular the energy and secondary particle phase space regions that are most important 
for the observed characteristics of EAS. 

The electromagnetic component of a shower is well determined by the depth of max¬ 
imum and the energy of the shower. Due to the electromagnetic cascade, having a short 
radiation length of ~ 36g/cm 2 , any information on the initial distribution of photons pro¬ 
duced in 7 r° decays is lost. Therefore the electromagnetic shower component depends on 
the primary particle type only through the depth of shower maximum. In contrast, the 
muon component is very sensitive to the characteristics of hadronic interactions. Once 
the hadronic shower particles have reached an energy at which charged pions and kaons 
decay, they produce muons which decouple from the shower cascade. The muons prop¬ 
agate to the detector with small energy loss and deflection and hence carry information 
on hadronic interactions in EAS. Due to the competition between interaction and decay, 
most of the muons are decay products of mesons that are produced in low-energy interac¬ 
tions. Therefore it is not surprising that muons in EAS are particularly sensitive to hadronic 
multiparticle production at low energy |:2][3). Recent model studies show that even at ultra- 
high shower energies the predictions on the lateral distribution of shower particles depend 
strongly on the applied low-energy interaction model mm 
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2 Muon production in extensive air showers 

Motivated by the measurement conditions of the KASCADE array 0, we consider 
showers with a primary energy of 10 15 eV and apply a muon detection threshold of 250 MeV. 
Using a modified version of the simulation package CORSIKA 0 we have simulated two 
samples of 1500 vertical and inclined (60°) proton and 500 iron induced showers. Below 
80 GeV the low-energy hadronic interaction model GHEISHA 2002 (8} and above 80 GeV 
the high-energy model QGSJET 01 (9) are applied. In the following vertical as well as 
60° inclined proton showers are discussed. The results are very similar for iron induced 
showers. 

In Fig. 0 the energy distribution of muons at detector level (1030g/cm 2 ) is shown for 
several lateral distance ranges. The maximum of this distribution shifts to lower energies 
for larger lateral distances. Most likely five consecutive hadronic interactions (number of 
generations) take place before a hadron decays into a muon, see Fig. [2] Here and in the 
following we consider only those muons that reach the ground level with an energy above 
the detection threshold. The number of generations show only a weak dependence on the 
lateral distance. 

To study the hadronic ancestors of muons in EAS, we introduce the terms grandmother 
and mother particle for each observed muon. The grandmother particle is the hadron in¬ 
ducing the last hadronic interaction that finally leads to a meson (mother particle) which 
decays into the corresponding muon. Most of the grandmother and mother particles are 
pions, but also about 20% of the grandmother particles are nucleons and a few are kaons. 
Details of the composition of mother and grandmother particles are given in Tab.0 



Fig. 1. Simulated energy distribution of muons 
for different lateral distances in vertical proton in¬ 
duced showers with a primary energy of 10 lo eV. 



Fig. 2. Number of generations before producing 
a muon visible at ground level (shown for various 
lateral distances). 


3 Relevant energy range 

The energy spectra of different grandmother particles produced in vertical proton show¬ 
ers are shown in Fig.0(left). They cover a large energy range up to the primary energy with 
a maximum at about 100 GeV. The maximum shifts to higher energies (several lOOGeV) 
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Table 1. Particle types of mother and grandmother particles in a vertical proton induced shower at 

10 15 eV. 



mother 

grandmother 

pions 

89.2% 

72.3% 

kaons 

10.5% 

6.5% 

nucleons 

- 

20.9% 


for 60° inclined proton showers, see Fig.^(left). The peak at 10 6 GeV in the nucleon en¬ 
ergy spectrum shows that also a fraction of muons stems from decays of mesons produced 
in the first interaction in a shower. Furthermore, the step at 80 GeV clearly indicates a mis¬ 
match between the predictions of the low-energy model GHEISHA and the high-energy 
model QGSJET. In Fig. 0 (right) the grandmother particle energy spectrum is shown for 
different ranges of lateral muon distance. The maximum shifts with larger lateral distance 
to lower energies. The same behaviour is visible for inclined showers, see Fig. 0] (right). 

Comparing the last interaction in EAS with collisions studied at accelerators, one has 
to keep in mind that the grandmother particle corresponds to the beam particle and the 
mother particle is equivalent to a secondary particle produced in e.g. a minimum bias p-N 
interaction. The most probable energy of the grandmother particle is within the range of 
beam energies of fixed target experiments e.g. at the PS and SPS accelerators at CERN. 




Fig. 3. Energy distribution of grandmother particles in vertical proton showers. Left panel: different 
grandmother particle types for a muon lateral distance range of 0-500 m at ground level. Right panel: 
different lateral distances, all particle types are summed up. 


4 Relevant phase space regions 

The further study of the relevant phase space of the mother particles is done for two 
different grandmother energy ranges and muon lateral distance ranges at ground level, see 
Tab. |3 The lateral distance ranges are chosen to resemble typical lateral distances measured 
at KASCADE and KASCADE-Grande, respectively ITOl . Motivated by the availability of 
protons as beam particles at accelerators we consider only those last interactions in EAS 
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Fig. 4. Same plots as in Fig|3] but for 60° inclined proton showers. 

that are initiated by nucleons. Exemplary, the results of this analysis are only shown for 
vertical proton showers. 

In Fig. [5] the rapidity spectra of mother particles (left: pions, right: kaons) are com¬ 
pared to the spectra of secondary particles of minimum bias proton-carbon and proton-air 
collisions with a fixed energy simulated with QGSJET labeled as fixed target. The spec¬ 
tra of mother particles in air showers are scaled to fit the falling tail of the fixed-energy 
collision spectra. No significant differences are found comparing the rapidity distributions 
of secondary particles in proton-carbon and proton-air collisions. As a consequence of the 
different selection criteria, the forward hemisphere in the mother rapidity spectra is clearly 
favoured compared to the spectra of secondaries in minimum bias collisions. The reason 
for this behaviour is the fact that the secondary particles (pions, kaons) are measured di¬ 
rectly in fixed target experiments whereas in EAS an additional condition for the mother 
mesons is applied. In order to get the information of the mother meson, it has first to decay 
into a muon which is detectable. At low rapidity the mother mesons are missing in the EAS 
distribution because the energy of the daughter muon is lower than the applied detection 
threshold or the muon decays and it is not visible. The reason for the missing pions in the 
EAS distribution at the highest rapidities is the interplay between decay and further inter¬ 
actions, that means the pions do not decay but undergo further interactions. This effect is 
not much pronounced for the kaon distribution, because the energy, for which the proba¬ 
bility of decay and interaction is of the same order, is eight times higher for kaons than for 
pions. 


Table 2. Energy and lateral distance ranges used for this analysis. 


energy range 

average energy 

lateral distance range 

80-400GeV 

160 GeV 

50-200 m 

30-60 GeV 

40 GeV 

200-600 m 
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Fig. 5. Rapidity distributions of mother particles (filled curves) compared with rapidity distributions 
of secondary particles in simulated single p+C (solid line) and simulated p+air (dashed line) colli¬ 
sions. Left panel: pions, right panel: kaons. The energy range of the grandmother particle is limited 
to 80-400 GeV and the lateral distance of the muons to 50-200 m to match experimentally accessible 
regions. The fixed target collision simulation is done at 160 GeV, corresponding approximately to the 
mean grandmother energy. The rapidity is normalized to the rapidity of the beam and grandmother 

particles, respectively. 

The phase space regions of mother particles produced in the last interaction in EAS are 
shown in Fig. ©and© We choose a nucleon as the grandmother particle with a mean energy 
of 160GeV (Fig.© and 40GeV (Fig. 13 respectively. The transverse momentum of the 
mother particle is plotted vs. its rapidity divided by the rapidity of the grandmother particle 
which is equivalent to the beam rapidity in fixed target reactions. On the left hand side this 
distribution is shown for pions as mother particles, on the right hand side for kaons. The 
maximum of these distributions, which shows the most important phase space region for 
the muon production in EAS, is at p± « 0.1 GeV and 0.7 in relative rapidity units (for 
a mean grandmother energy of 160 GeV) and shifts to slightly higher p± but stays at the 
same rapidity for a mean grandmother energy of 40 GeV. In both cases the distributions 
of pions and kaons are similar. For kaons higher particle transverse momenta are more 
important than for pions. The phase space regions of relevance to EAS are summarized in 
Tab.©and indicated with the dashed (red) boxes in Fig. ©and© 

Table 3. Phase space regions of hadronic interactions relevant for muon production in EAS. 


average energy (GeV) 

V/ybeam 

p± (GeV/c) 

160 

0.3 - 1.1 

0.0 - 0.7 

40 

0.3 - 1.1 

O 

l 

o 

o 
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Fig. 6. Phase space of mother particles. Grandmother energy range: 80-400 GeV. Left panel: pions, 
right panel: kaons. The dashed box (red) indicates the most interesting phase space region which 
includes more than 90% of this particles. 




Fig. 7. Same plots as in Fig. [6] but the grandmother energy range is 30-60 GeV. 


5 Existing fixed target data 

In contrast to collider experiments, in fixed target setups it is more straight forward 
to measure particles in forward direction. The energy range of most of the existing fixed 
target experiments has a broad overlap with the energy range of relevance to EAS. 

There exist a number of fixed target measurements employing a proton beam with low 
energy (<20GeV) with a good phase space coverage, however, data at higher energy and 
for pion projectiles are very sparse. 

In Fig. 0a compilation of existing fixed target data is shown. Mostly the data on sec¬ 
ondary particle distributions are given in the observables angle between beam and sec¬ 
ondary particle momentum (0) and the momentum of the secondary particle itself (p sec ). 
Some measurements are also published using rapidity (y) and transverse momentum of 
secondary particles (p±). For an easier comparison, in Fig. 0all data sets are shown in the 
rapidity-transverse momentum plane, whereas the phase space regions covered by 0 — p sec 
measurements have been converted approximately. 
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In both figures the existing data are indicated by shaded (colored) regions, whereas the 
beam momentum is given as abscissa and the secondary particle observable as ordinate. 
The relevant phase space is shown as box histogram. In this case the grandmother momen¬ 
tum is used instead of the beam momentum for the abzissa. The size of the boxes indicates 
the relative importance of the beam and secondary particle momenta. 

The work of Eichten et al. CD has become a widely used reference data set. This 
experiment used a proton beam with a beam momentum of 24 GeV and a beryllium tar¬ 
get. The secondary particles (pions, kaons, protons) are measured in a broad angle range 
(17mrad < 9 < 127mrad) and in a momentum region from 4GeV/c up to 18GeV/c. 
Other measurements cover only a smaller part of the phase space of interest to EAS 
1121 1131 ITll n~5l Uhl ITT! II81 1191 . A measurement for preparing the CERN neutrino beam 
experiments was performed by SPY/NA56 (NA52) [20, 2JJ. Therefore a 450 GeV proton 
beam and a thick beryllium target was used. Because of the very limited angular range 
(9 « 0°) it is not included in Fig.[8]or Fig.[9] 
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Fig. 8. Coverage of the phase space regions of relevance to EAS (box histograms) by existing 
fixed target data using a proton beam and a beryllium or carbon target (shaded/colored regions). Left 
panel: total momentum of secondary pions vs. total momentum of proton projectiles. Right panel: 
angle between beam and secondary particle momentum vs. beam momentum. 


It is clear from this comparison (Figs.[8]an d0 that existing data do not cover the most 
important energy and secondary particle phase space regions. However, there are three 
promising current ongoing data analyses. The HARP (PS214) experiment t22l at the PS 
accelerator at CERN has taken p+C and 7T+C data at 15 GeV in 2001. The MIPP (FNAL- 
E907) experiment 1231 at Fermilab just started data taking in the beam energy range 5- 
120GeV in 2005. Because of its excellent forward acceptance, it is also thought about 
using the NA49 detector l24l at the SPS accelerator at CERN in a follow-up experiment to 
perform similar measurements at higher energies (up to 160 GeV) G3 
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Fig. 9. Compilation of the phase space regions covered by fixed target data given in transverse 
momentum and rapidity of secondary particles and the phase space regions covered by the 8 — p sec 
data (see Fig.[8), whereas an approximate conversion of the covered phase space has been done. Left 
panel: transverse momentum of secondary pions vs. total momentum of proton projectiles. Right 
panel: rapidity of secondary pions normalized by the beam rapidity vs. beam momentum. 


6 Conclusions and outlook 


Due to the interplay between decay and interaction of pions and kaons, low energy 
hadronic interactions are very important for muon production in extensive air showers. 
With increasing lateral distance the mean energy of these interactions, which are mainly 
initiated by pions and nucleons, decreases. 

The most important interaction energies and phase space regions are accessible for 
fixed target experiments. However, so far only measurements with protons as projectiles 
and with a very limited secondary particle phase space exist. The situation could be im¬ 
proved considerably, if data from fixed target experiments with large acceptance detectors 
such as HARP, NA49 and MIPP would be analysed for minimum bias collisions especially 
for 7T+C reactions. 

Fixed target measurements can contribute to improving low energy interaction models. 
Due to the sensitivity of EAS to low energy hadron production, such improvements will 
increase the reliability of air shower simulations. 

Finally it should be noted that multi-detector installations such as KASCADE can be 
used to check the consistency of the simulation of hadronic interactions by comparing 
different observables of EAS and their correlation with model predictions. 

Acknowledgements: The authors thank Dieter Heck for many fruitful discussions and 
help with modifying CORSIKA to include the muon ancestor information. 
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